Is cosmic acceleration slowing down? 
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We investigate the course of cosmic expansion in its recent past using the Constitution SN la 
sample (which includes CfA data at low redshifts), jointly with signatures of baryon acoustic os- 
cillations (BAO) in the galaxy distribution and fluctuations in the cosmic microwave background 
(CMB). Earlier SN la data sets could not address this issue because of a paucity of data at low 
redshifts. Allowing the equation of state of dark energy (DE) to vary, we find that a coasting model 
of the universe [qo = 0) fits the data about as well as ACDM. This effect, which is most clearly seen 
using the recently introduced Om diagnostic, corresponds to an increase of Om and q at redshifts 
z < 0.3. In geometrical terms, this suggests that cosmic acceleration may have already peaked 
and that we are currently witnessing its slowing down. The case for evolving DE strengthens if 
a subsample of the Constitution set consisting of SNLS+ESSENCE+CfA SN la data is analysed 
in combination with BAO+CMB using the same statistical methods. The effect we observe could 
correspond to DE decaying into dark matter (or something else). A toy model which mimics this 
process agrees well with the combined SN la+BAO+CMB data. 

PACS numbers: 98.80.Es,98.65.Dx,98.62.Sb 



The existence of cosmic acceleration at redshifts less 
than ~ 0.5 appears to be well established by several inde- 
pendent data sets including: SN la luminosity distances, 
cosmic microwave background temperature anisotropy 
and polarization maps, and baryon acoustic oscillations 
in the galaxy power spectrum. Most recent analysis 
performed using the data mentioned above [ij, as well 
as data from Chandra and SDSS 0] cluster cata- 
logues, show that if the DE equation of state (EOS) 
w = pde/pde is assumed to be a constant, then there 
remains little room for departure of DE from the cosmo- 
logical constant, since < 0.06 at the la confidence 

level. However, in the absence of compelling theoretical 
models with an unevolving EOS, one must reexamine the 
data impartially by removing this prior, if one is to look 
for serious alternatives to the cosmological constant 

A large sample of nearby SN la with z < 0.08 has 
recently been published [^. Adding this to the Union 
sample @, leads to the Constitution set Q which is the 
largest SN la sample to date. Consequently one might 
expect significantly better limits on 1 -t- w to be placed 
from this new sample. As shown in Q the Constitu- 
tion set with a BAO prior leads to approximately the 
same upper limit on |1 + wj as that obtained from all 
other data sets if u; = const is assumed. However, as 
noted earlier, the assumption of a constant EOS is not 
very realistic. For this reason, in this letter, we drop this 
assumption when analyzing data from the Constitution 
set, together with BAO data at z = 0.2 and z = 0.35 @ 
and the observed CMB shift parameter R obtained from 
acoustic oscillations in the CMB temperature anisotropy 
power spectrum ll| . The data is analyzed using the pop- 
ular CPL ansatz together with the recently introduced 




FIG. 1: The Om diagnostic is shown as a function of 
redshift for DE models with Qom ~ 0.27 and w = 
-1,-0.8,-0.6,-0.4,-0.2 (bottom to top). For Phantom 
models (not shown) Om would have the opposite curvature. 
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where 



flarnil + zf + f^DE 



n^E = (l-f7o™)exp|3 j^'ii^^dzj(2) 

is the expansion history of a spatiaUy flat Friedmann- 
Robertson- Walker (FRW) Universe with scale factor a{t) 
and Hubble parameter H(z) = a/a. The CPL ansatz ex- 
presses the EOS in terms of the redshift z in the following 
form 0: 



w{z) — Wq + Wi 
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In contrast to w{z) and the deceleration parameter 
q{z) = —d/aH^, the Om{z) diagnostic depends upon 
no higher derivative of the luminosity distance than the 
first one. Therefore, it is less sensitive to observational 
errors than either w or q. Om is also distinguished by the 
fact that Om{z) = ilom for ACDM . Om is very useful 
in establishing the properties of DE. For an unevolving 
EOS: l + w~ [Om(z) - f^om](l - f^om)^^ at z < 1, 
consequently a larger Om{z) is indicative of a larger w; 
while at high z, Om{z) fiom, as shown in figure [TJ 

The present analysis uses the recently compiled 'Con- 
stitution set' of 397 type la supcrnovae covering a red- 
shift range from z^in = 0.015 to z^^ax = 1.551. The Con- 
stitution set is the largest SN la luminosity distance sam- 
ple currently available and includes 139 SN la at 2 < 0.08. 
Our analysis considers the SN la data individually as well 
as in combination with BAO distance measurements ob- 
tained at z — 0.2 and 2 — 0.35 from the joint analysis of 
the 2dFGRS and SDSS data 3. The BAO distance ratio 
Dv(z = 0.3b) /Dv{z = 0.20) = 1.736 ± 0.065 was shown 
in [8| to be a relatively model independent quantity. Here 
Dv{z) is defined as 
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We also use the CMB shift parameter [1] which is the re- 
duced distance to the last scattering surface (zjg = 1090) 



dz 



1.71 ±0.019 



(5) 



While SN la and BAO data contain information about 
the Universe at relatively low redshifts, the R parameter 
probes the entire expansion history up to last scattering 
at zis. Because of this, our analysis also examines the 
goodness of fit for CPL parametrization when applied 
simultaneously to data at low and high redshifts. 

Fig. [2] shows q{z) and Om{z) reconstructed using ([3]) 
and SN la -I- BAO data (upper panels) and SN la -I- 
BAO CMB data (lower panels). Red lines are best fit 
reconstructions while dashed green lines show la CL's 



from a analysis. It is interesting that the best fit fiat 
ACDM model (fJo„ = 0.287) satisfying SN la-^BAO data 
does not lie within the la CL of our best reconstruction. 
By contrast, fT] obtain l+w = O.OlStpQgg after assuming 
w — constant, which underscores the difference made by 
dropping the w — constant constraint. It is interesting to 
mention that the reduced also drops from x^^^ = 1.182 
in case of ACDM model to Xred = 1-171 in case of varying 
dark energy model which makes the assumption of the 
additional parameter worthwhile. 

The growth in the value of Om{z) at low z in the upper 
panel of fig.[2]is striking, and appears to favor a DE model 
with an EOS which increases at late times. This could be 
preliminary evidence for a decaying DE model since, in 
this case, the EOS would increase at late times, resulting 
in an increase in the low 2 value of the Hubble parameter 
and therefore also of 0771(2). 

These results change dramatically with the inclusion 
of CMB data. The lower panel of fig. [2] shows that our 
reconstruction of Om{z) is now perfectly consistent with 
ACDM - for which Om{z) is unevolving. However, one 
should not stop and conclude that the standard model is 
confirmed once more. Rather, what one observes here is 
an incompatibility of the CPL parametrization for ■w{z) 
with this combination of data sets. In other words, the 
functional form of the CPL ansatz is unable to fit the 
data simultaneously at low and high redshifts. This can 
be clearly seen if we compare Xsn+bao ~ 461.63 {xled ~ 
1.171) for the best fit obtained using SN+BAO data, 
with the significantly larger x'sn+bao+cmb = 

467.07 

{xled = 1-182) obtained using SN+BA0H-CMB(R) data. 
The addition of one more data point (the CMB shift pa- 
rameter), increases the best fit x'^ by niore than 5. Sup- 
port for this viewpoint is also provided by fig. [3l which 
shows the best fit regions in parameter space obtained us- 
ing the CPL ansatz after fitting to SN la (red pluses), SN 
la+BAO (green crosses) and SN la-HBAO-f CMB (blue 
stars) data. Contrast the good overlap between best fit 
regions obtained using SN la and SN Ia-|-BAO data, with 
the relative isolation of the best fit region obtained using 
SN la-HBAO-t-CMB data. 

Our reconstruction of Om{z) appears to favor DE with 
an increasing EOS at low redshifts 2 <, 0.3 (fig. [2]upper- 
right panel). We have also seen that the CPL ansatz is 
strained to describe the DE behavior suggested by data 
at low and high z. We believe the reason for this stems 
from the fact that the CPL ansatz implicitly assumes that 
the redshift interval from 2 = to 2 ~ 2 represents noth- 
ing special for DE, so that w{z) can safely be expanded 
in a Taylor series in powers of 2/(1 -I- 2) in this inter- 
val. This need not, however, be true, since models have 
been suggested in which dark energy decays with a char- 
acteristic time of order of the present age of the Universe 
0, [13, ■ Note also that the large negative value of wi 
(fig. [3]) suggests that DE was practically nonexistant at 
high redshifts. This too could be an ansatz-related fea- 
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FIG. 2: Reconstructed q{z) and Om{z) from SN la + BAO data (upper panels) and SN la + BAO + CMB data (lower panels) 
using the GPL ansatz. Solid red lines show the best fit values of Om{z) and q{z) while dashed green lines show the la GL. 
The dramatic difference between the upper panels and the lower one is indicative of the inability of the GPL parametrization 
to fit the data at low and high redshifts simultaneously. The spatially flat AGDM model corresponds to a horizontal line with 
Om{z) = ilom in the right panels (not shown). 



ture, since the reliability of ([3]) at high redshifts remains 
somewhat ambiguous. Keeping these issues in mind, we 
re-analyze the data using a simple toy model for w{z) 
which encapsulates the main features of the effect discov- 
ered, and show that this ansatz can provide a better fit 
to the combination of SN-hBAO-hCMB data than CPL. 
Our ansatz is 



w{z) = — 



l + tanh[(z-Zt)A] 



(6) 



(a similar form was used in 1J| and other papers to search 
for fast phase transitions in DE at larger values of z). 
This fit ensures w = —1 at early times, and then increases 
the EOS to a maximum of w ~ at low z. It has the same 
number of free parameters as the CPL ansatz but does 
not permit the crossing of the phantom divide at w = —1. 
The best fit cosmology obtained using this ansatz has 



zt = 0.008 {i.e. zt « 0), A = 12.8, f7o,„ = 0.255 and 
XsN+BAO+CMB ^ 466.50, and presents an improvement 
(Ax^ = —0.6) over the best fit for the same data set ob- 
tained using the CPL ansatz. We would like to emphasize 
that this new fit is only used as an example to demon- 
strate that the CPL ansatz may not be flexible enough to 
determine cosmological parameters for a class of rapidly 
evolving DE models. Fig. |4] shows the deceleration pa- 
rameter q and the Om diagnostic reconstructed using Q 
& Interestingly, ACDM as well as a universe which is 
currently coasting [qo ~ 0), can both be accommodated 
by the data at roughly the same level of confidence ! 

To check the robustness of our results we redo our anal- 
ysis on a subsample of the Constitution data set, namely 
the SNLS+ESSENCE+CfA SN la data (234 data points 
in all). Our results, shown in figurelHlfor the CPL ansatz, 




FIG. 3: la contours for GPL parameters uiq-uii (left panel) and luo-i^om (right panel) reconstructed using SN la data (red 
pluses), SN la+BAO data (green crosses) and SN la+BAO+GMB data (blue stars). Note the consistency between SN la 
and BAG data. The absence of any overlap between the la contours for SN la+BAO+GMB and SN la+BAG data could be 
indicative of tension between the GPL parametrization and the data. We should note here that plotting small contours that 
indicate tight constraints on the parameters is misleading when the best fit result produces a bad fit to the data (blue contours) . 
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FIG. 4: The cosmological deceleration parameter q{z) (left panel) and Om(z) (right panel) reconstructed using a combination 
of SN la, BAG and GMB data and the ansatz ((6|. Solid red lines show best fit reconstructed results while dashed green lines 
show reconstructed results within la GL. 



are summarized below: 

1. An excellent overlap exists between the la con- 
tours wq-wi (top left) and WQ-ilom (top right) re- 
constructed using SN la, SN la-hBAO and SN 
la+BAO+CMB data. This demonstrates that the 
CPL ansatz works quite well for this combination 
of data sets and hints that the tension noticed 
in figure [3] could be coming from data sets which 
have been excluded from the present SN la com- 
pilation: namely the Gold data, the high z HST 



data and older SN la data sets. The visual im- 
pression conveyed by this panel, which appears to 
support evolving DE, receives statistical support: 
the best fit for SN la data is 267.69, while 
= 267.92 for SN Ia-|-BAO and = 268.89 for 
SN la-hBAO-l-CMB. We therefore find that the x^ 
values for the three data sets (SN la, SN la+BAO, 
SN la-l-BAO-l-CMB) he much closer together for 
the data shown in figure El compared to the data 
in figure [31 




FIG. 5: Top: la contours for GPL parameters wq-wi (left panel) and wo-Qom (right panel) reconstructed using SN la data 
(SNLS+ESSENCE+GfA, red pluses), SN la+BAO data (green crosses) and SN la+BAO+GMB data (blue stars). Note the 
compatibility between the different data sets and that spatially flat AGDM (red cross at wo = — = in the left panel) 
appears to be in tension with this combination of data. Middle and Bottom: Reconstructed q{z) and Om{z) from SN la + 
BAG data (middle). Reconstructed q{z) and Om{z) from SN la + BAG + GMB data (bottom). In the middle and bottom 
panels solid red lines show the best fit values of Om{z) and q{z) while dashed green lines show the la GL. In the bottom panel, 
the best fit values for AGDM are: qo ~ —0.6 and Om(z) — 0.255. All results have been obtained using the GPL ansatz ([Sj. 



2. A larger value of Om{z) and q{z) at low red- 



shifts is supported by the present analysis of 
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SNLS+ESSENCE+CfA supernovae in combina- 
tion with BAO and CMB data (figure [5] middle 
and bottom). Indeed, coasting cosmology (go — 0) 
provides an excellent fit to the data while ACDM 
(qo ~ —0.6) appears to be excluded at la by this 
data set. This is in marked contrast to the results 
in fig. [2] (lower panel) which distinctly favoured 
ACDM . 

3. The fact that the spatially flat ACDM shows 
weaker consistency with the SN la subsample + 
BAO + CMB data is clearly seen from the best 
fit values for ACDM : (i) = 274.64 using only 
SN la {flm = 0.28), (ii) = 275.87 for SN 
Ia+BA0(f7,„ = 0.28), (iii) = 276.84 for SN 
la+BAO+CMB {n„, = 0.26). Comparing with 
the results for evolving DE discussed earlier, we 
find that the incremental value of Ax^ between 
the best fit evolving DE model and best fit ACDM 
is ~ 8, and favours evolving DE (the reduced x^ 
drops from xled — 1-188 in case of ACDM model to 
Xred ~ 1-159 in case of varying dark energy model). 

To summarize, the recently released Constitution SN 
la data set appears to support DE evolution at low red- 
shifts. There also appears to be some tension between 
low z (Constitution SN la + BAO) and high z (CMB) 
data, when analyzed using the CPL ansatz. (However, 
this tension decreases when only a subsample of the Con- 
stitution set is analysed.) There could be several reasons 
for this. 

• Systematics in some of the data sets is not suffi- 
ciently well understood. This effect may have a 
purely astronomical explanation and be a result 
of some systematic effect, e.g. if new nearby CfA 
SN la are brighter on average. However, we have 
found additionally found that if the effect is as- 
sumed to be cosmological, then the implied DE 
behavior at low redshifts (using SN la data) is 
more consistent with a rather large value of the 
Dviz = 0.35)/ Dv{z = 0.20) BAO distance ratio 
derived in [§]. 

• Another possibility is that this behavior of 
Om, w, q is an apparent one, which is induced by 
a local spatial inhomogeneity - a kind of a "Hubble 
bubble" , but with a large-scale matter overdensity 
(though the results of 7\ do not appear to support 
its existence). 

• Different SN la subsamples comprising the Consti- 
tution set have varying properties. 



suggested by the data. 
Clearly one must wait for more data before deciding 
between these alternatives. Note that the bulwark of sup- 
port for evolving DE comes from BAO data in conjunc- 
tion with Constitution SN la data [3] which includes 139 
SN la at z < 0.08. It is therefore encouraging that several 
hundred SN la light curves are expected to be presented 
by the Nearby Supernova Factory (0.03 < z < 0.08) and 
SDSS (0.05 < z < 0.35). As concerns ACDM, our con- 
clusion is that it is still viable at the 95% confidence level 
but perhaps one should take its subtle challenges more 
seriously [l5|. 
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• The CPL ansatz is not versatile enough to accom- 
modate the cosmological evolution of dark energy 



